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Theoretical studies on the possible origin of room temperature ferromagnetism (ferromagnetic 
once crystallized) in the molecular transition metal complex, V(TCNE) :E (1.5 < x < 2) have 
been carried out. For this family, there have been no definite understanding of crystal structure 
so far because of sample quality, though the effective valence of V is known to be close to +2. 
Proposing a new crystal structure for the stoichiometric case of x=2, where the valence of each 
TCNE molecule is -1 and resistivity shows insulating behavior, exchange interaction among d- 
electrons on adjacent V atoms has been estimated based on the cluster with 3 vanadium atoms 
and one TCNE molecule. It turns out that Hund's coupling among d orbitals within the same 
V atoms and antiferromagnetic coupling between d oribitals and LUMO of TCNE (bridging V 
atoms) due to hybridization result in overall ferromagnetism (to be precise, ferrimagnetism). 
This view based on localized electrons is supplemented by the band picture, which indicates the 
existence of a flat band expected to lead to ferromagnetism as well consistent with the localized 
view. The off-stoichiometric cases [x < 2), which still show ferromagnetism but semiconducting 
transport properties, have been analyzed as due to Anderson localization. 
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1. Introduction 

Recently, transition metal complexes with organic 
molecules have attracted much attention. In particular, 
" ferromagnetism" was found to be realized in some com- 
pounds, such as T„(C6H 6 )„ + i (T=Sc, Ti, and V) 1 ) and 
V(TCNE) ;C (TCNE= tetracyanoethylene, 1.5 < x < 2) 
.2-4) ^ ne i a ^ er compound has been discussed intensively 
for the past two decades since its discovery in 1991, 2 ) be- 
cause of its high Curie temperature T c ~ 400K, the high- 
est of the transition-metal-organic-molecule complexes. 
However, the mechanism of its ferromagnetism has not 
been clarified yet. This is partly because the crystal 
structure of these materials has not been identified so 
far, while many crystal structures of V(TCNE)a; have 
been proposed. On the basis of the band structure calcu- 
lations, the stability of ferromagnetism was discussed. 5 ' 6 - ) 
Nevertheless, the essence of the mechanism of ferromag- 
netism is not yet fully understood, because these calcu- 
lations treat many body effects only in onc-clcctron ap- 
proximation and then effects of strong correlations have 
not been fully assessed. 

The purpose of this paper is twofold. First, we pro- 
pose a new structure of V(TCNE) :E which satisfies all 
the experimental constraints. Second, on the basis of this 
structure, 1) the exchange interaction among electrons 
at adjacent V ions is estimated by numerical diagonal- 
ization method which can take the electron correlation 
fully into account, and 2) the effect of the itinerancy of 
electrons is discussed by the tight-binding model, which 
shows that a sharp peak of the density of states (DOS) 
at the Fermi level appears due to the existence of the flat 
band when the couplings between two dimensional layers 



are absent. The latter result can well explain the trans- 
port properties of the system exhibiting a resistivity of 
the variable-range-hopping behavior. 

2. Electronic State of TCNE 

It is reported that the lowest unoccupied molecular 
orbital (LUMO) of TCNE strongly hybridizes with d or- 
bitals of V. 7 ' The electronic state of LUMO has been 
discussed on the basis of the quantum chemical calcula- 
tion. 8 ) In this section, we show that LUMO is expressed 
qualitatively in the 7r electron approximation. 

The structure of TCNE is shown in Fig. 1(a) where 
Cj and Nj are i-th Carbon (C) and Nitrogen (N). Three 
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Fig. 1. (a) Structure of TCNE. The distance between Ns is about 
AAA. (b) Energy level scheme of the molecular orbitals of TCNE 
on the 7r electron approximation. 



transfer integrals t, t' and t" among 7r-orbitals of C and 
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N are shown in Fig. 1(a). Fig. 1(b) shows the energy level 
scheme of the molecular orbitals of TCNE in the tight- 
binding approximation for the resonable values of the 
transfer parameters t' = l.lt and t" = 1.5i, because the 
transfer integrals are propotional to d~ 2 where d is the 
distance between C-C (d ~ 1.431), C=C (d ~ 1.371) 
and C=N (d ~ 1.161). 8 > Then, the wave function of 
LUMO $l is written as $j, = Ai(—ip Nl — ip N2 + ip N3 + 

VnJ + A 2 {tp Cl + fc 2 ~ Vcs ~ <Pc 6 ) + A 3 ((pc 3 - fc A ) 
where ip^i and ifd are the wave functions of 7r-orbital of 
i-th N and C, respectively, and A x = 0.355, A 2 = 0.129 
and A3 = 0.463 are obtained with the same parameter 
set. LUMO is consistent with the wave function on the 
basis of the quantum chemical calculation. 8 ) LUMO has 
a large weight of it orbitals of N. Namely, it suggests the 
imporancc of the hybridization between tt orbitals of N 
and d orbitals of V. 

It is also noted that the transfer integrals are gener- 
ally propotional to d~ 4 . 9 ^ At that time, these transfer 
integrals and the HOMO-LUMO gap are estimated as 
t' = 1.2t and t" = 2.3t, and of the order of It. The weight 
of C and N in LUMO is a little different from the case 
that the transfar integrals are propotinal to dr 2 . How- 
ever, the energy level scheme of the molecular orbitals is 
qualitatively consistent with the case of d~ 2 . 

3. Structure of V(TCNE) X 

We propose a new structure of V(TCNE).e on the ba- 
sis of experimental results. The coordination number of 
V in V(TCNE) :E -yS and V(TCNE) X has not been de- 
termined, because of the lack of x-ray diffraction lines 
due to the existence of defects of V. 4 ) Nevertheless, the 
number of ligands around V has been studied in detail 
by EXAFS, 10 ) showing that V has a slightly distorted 
octahedral environment constructed by six N atoms. It 
was also proposed that these N atoms belong to differ- 
ent TCNEs from the measurement of the N-N distance 
of octahedra (2.8 1) which is shorter than the distance 
between two Ns in the TCNE molecule (4.4 1). 

The nominal valence ny of V is estimated as ny = +2 
by XANES. 10) At that time, the n v /6 electrons of V 
transfer to one TCNE, because V bonds with six TCNEs. 
The charge neutrality requires that the valence utc ne of 
TCNE should be Utcne = —ny/x. Then, the number 
Ny of Vs bonding with TCNE is given as Nyny/6 = 
ny/x, and corresponds to Ny = 3 for x = 2 and Ny = 4 
for x = 1.5. 

Based on these experimental results, we propose the 
structure of V(TCNE) :E as shown in Fig. 2: Fig. 2(a) 
for V(TCNE)i. 5 and Fig. 2(c) for V(TCNE) 2 . While 
the structure of Fig. 2(a) has already been proposed, 4 ) 
the structure of V(TCNE)2 is new. The structure of 
V(TCNE) X , for 1.5 < x < 2, is expected to be an in- 
termediate one between Fig. 2(a) and Fig. 2(c), which is 
obtained by extracting V ions from Fig. 2(a), as shown 
in Fig. 2(b). Note that Vs are randomly extracted from 
V(TCNE)i.5, but all Vs have six-coordination. 

4. Local Electronic State of V 

To construct an effective model of V(TCNE) :E , it is cru- 
cial to understand the local electronic state of V. Since 



b) 1.5 < % < 2 




Fig. 2. Structures of V(TCNE) I . Fig. 1(a), (b) and (c) are the 
structure of x = 1.5, 1.5 < x < 2 , and x = 2, respectively. Red, 
blown and blue circles indicate V, C and N, respectively. 



V is surrounded by six N atoms forming a slightly dis- 
torted octahedral structure, we assume that the symme- 
try around V is D 4 /j, although the exact location of the 
apical N atoms has not been determined yet. The crys- 
talline electric field (CEF) levels of the d orbitals of V 
depends on the location of apical N atoms. 

From the symmetry requirement, the hybridizations 
are allowed only between d yz and d zx orbitals of V and 
7r orbitals of N. Here, we have defined the x and y axes so 
that N atoms of the V-TCNE plane will lie on the x and 
y axes, and V atoms at the origin. We also have defined 
the z axis so that the apical N atom will lie on the z axis. 
On the other hand, the d xy orbital of V can be regarded 
as an almost localized one. We also find that e g orbitals 
of V do not hybridize with n orbitals, but hybridize with 
a orbitals of N. The a bonding formed between these or- 
bitals makes a strong covalent bonding. 7 ) Namely, these 
bonding orbitals are important for stabilizing the struc- 
ture of VCTCNE)*. 

5. Cluster Model - Estimate of Exchange Inter- 
cation 

Here, we construct an effective cluster model for x — 2 
focusing on three Vs and TCNE as shown in Fig. 3(a). 
We also discuss the case of x=1.5 for four Vs and TCNE 
as shown in Fig. 3(b). We will estimate the exchange 
interaction among electrons at two adjacent V atoms. 
The Hamiltonian for the case of Fig. 3(a) consists of three 



a) V(TCNE) 2 b) V(TCNE) ]5 




Fig. 3. (a) Cluster of three Vs and TCNE. (b) Cluster of four Vs 
and TCNE. 



parts, the hybridization part H$, the d-site part and 
the LUMO part -Hlumo as 

H = H + Hd + i?LUMO- (1) 
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Explicit forms of these terms are written as follows: 
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exchange correlation. Similarly, spin correlation between 
electrons in d(2 tZX ) (or rf(3 i2X )) and LUMO becomes an- 
tiferromagnetic, causing a ferromagnetic correlation be- 
tween e?(i.j,z) and d(2, zx ) ( or ^(3,zs)) orbitals. Therefore, 
electrons hopping among d y2 (or d zx ) orbitals and LUMO 
takes spin state, S = 1 in the truncated system 
and 5* = in the full lattice systems of d yz and d zx 
orbitals and LUMO of V(TCNE) 2 . On the other hand, 
spins of d/i tVZ \, rf(i )ZX ) and d(i <xy > become parallel due 
to the Hund's rule coupling as shown in Fig. 4(a). Then, 
considering t ~ leV, the Curie temperature T c , which 
will be of the order of effective exchange coupling, J, is 
estimated to be of the order of 1000(K), which is consis- 
tent with experimental values of T c . 
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Here, d(ij)a is an annihilation operator of an electron 
with spin a on the d orbital j at site i ;d(j x)<d ^(i,2)a and 
d(j ,3) CT represent d zx , d yz and d xy , respectively ,while (/Jo- 
is that on LUMO; t is the transfer integral between d 
orbitals and LUMO; Ad, Urf, U^ and Jd are one-body 
level of d xy orbital, intra-, inter-, and exchange Coulomb 
interaction parameters; Al and Ul are one-body level of 
LUMO and Coulomb interaction parameter of LUMO. 

In the case of V(TCNE)2, each V has three electrons 
on d orbitals and TCNE has one electron on LUMO. 
Thus, the total number of electrons is 10 (= 3 x 3 + 1) in 
the truncated model. We have studied the ground state 
of this model of the cluster of V(TCNE)2 on the basis 
of the numerical diagonalization method. It turns out 
that the ground state of the cluster of V(TCNE)2 has 
S=4 state as shown in Fig. 4(a) for the parameter set 
U d /t = 4.0, U'Jt = 3.0, J d /t = 0.5, U L /t = 2, A d /t = 0.0 
and Al/1 — 2.0. Note that the ground state of S — 
4 prevails in a wide region of the parameter space of 
A^/t and A d /t although it is not shown here. For the 
parameter set A d /t = and A L /t = 2.0, the effective 
exchange coupling J between d orbitals is estimated as 
23 / 1 = E s=3 /t - Es=i/t ~ 0.08, where E s=3 and E s=i 
are the energies of the first excited state and the ground 
state. 

The nature of the spin state of d(i, yz ) orbitals of V 
and LUMO of TCNE can be seen from the spin— spin 
correlations {s yz s z v ) between d( liI/z ) and tp, and {s yz s z zx ) 
between dn yz ) and ^(i,zx) or dn, xv )- They are estimated 
as (s yz s^,) ~ — 0.0672; and (s yz s z x ) ~ 0.18i for the same 
parameter set. These results imply that the antiferro- 
magnetic correlation works between dn yz -) and (p, and 
the ferromagnetic correlation works between dn tyz -\ and 
d(x,zx) or d(i, xy ). 

The origin of the ferromagnetic exchange between d 
orbitals at adjacent sites is understood on the basis of 
the valence bond picture as follows: due to the kinetic 
energy gain, spin correlation between electrons in d(i,y 2 ) 
and LUMO becomes antiferromagnetic as in the super- 
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Fig. 4. Snapshot of the spin states of (a) V(TCNE) 2 , and (b) 
V(TCNE)i.5 on the valence bond picture. 



This mechanism of the ferromagnetic correlation on 
the valence bond picture, in which electrons in LUMOs 
have opposite spin polarization of d yz and d zx orbitals 
by the 7r-d hybridization, corresponds to the mechanism 
which Kanamori and Terakura discussed on the band pic- 
ture. 11 - 1 In other words, the magnetically ordered state 
can be regarded as a kind of ferrimagnetism as shown in 
Fig. 4(a). Such a ferrimagnetic nature has been pointed 
out by L. Valade and her coworkers on the basis of mea- 
surement of saturated magnetic moment. 12 ^ 

The origin of this ferromagentism may be understood 
based on the molecular orbital picture as follows. The 
molecular orbitals $i ~ $ 4 of the cluster of V(TCNE) 2 , 
which are the eigenstates of the Hamiltonan H eq. (2), 
are given, in order of increasing eigen values, by 



$ 2 = 



a(d( ltVZ ) + d(2,zx) + d( 3 , zx )) + fitp (5) 
1 



V2 



(d(2,zx) — d(s iZX )) 



$ 3 = ^(<W) 



l (3,zx)J 



(6) 
(7) 



a ' {d(l,yz) + d(2,zx) + <^(3,z:r) )+/?V, (8) 

where a, j3, a' and j3' are the normalization factors, and 
a a' and /3 ^ /?'. The molecular orbitals $i and $4 
are the bonding and antibonding molecular orbitals. The 
orbitals $2 and $3 are non-bonding molecular orbitals 
(NBMOs) consisting of only d components. Without in- 
teractions, $1 and $2 or $3 ($2 and $3 are degenerate) 
are doubly occupied. In the presence of interactions, eq. 
(3) and (4), the occupation number of $1 ~ $4 orbitals is 
altered so as to minimize the total energy. Fig. 5 shows 
the result of the number of electrons in the molecular 
orbitals $1 ~ $4 as a function of A^/t in the Hilbert 
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space of S z = 4 on the basis of the numerical diagonaliza- 
tion method with the same parameter set as Ud/t = 4.0, 
U'Jt = 3.0, J d /t 
A L /t = 2.0. Here 
the up spin (down spin) electrons of the m-th molecular 
orbital. The number of electrons in NBMOs (m = 2 and 



= 0.5, U L /t = 2, A d /t = 0.0 and 
(n m i) indicates the number of 




Fig. 5. Number of electrons of the up spin (black line) or down 
spin (red line) in MOs of the cluster of V(TCNE)2 with the set of 
the parameters, U d /t = 4.0, U'Jt = 3.0, J d jt = 0.5, A d /t = 0.0, 
and Uj^/t = 2.0. Index m is the number of the molecular orbitals. 



3) is almost equal to one. Furthermore, the electrons in 
NBMOs are in a nearly full spin polarized state, i.e. "fer- 
romagnetic state" . On the other hand, the total numbers 
of electrons in $i and $4 are nearly two and S z ~ 0. The 
origin of the ferromagnetic correlation is also understood 
by the same way of the valence bond picture. 

In the case of V(TCNE)i.s, in which the total number 
of electrons is 13 (=3x4 + 1), the ground state has S = 
11/2 state in the truncated model such as V(TCNE)2- 
These results for the cluster of V(TCNE)i. 5 show that 
the origin of ferromagnetic correlations of V(TCNE)i.5 
is the same as in V(TCNE) 2 . 

Precisely speaking, the mechanism of the ferromagen- 
tism discussed here can be applied only to the insu- 
lating state where the electrons are localized on each 
orbital. However, V(TCNE) X (1.5 < x < 2) exhibits 
a temperature dependence on the conductivity a oc 
T- 1 / 2 exp[-(To/T) 1 / 4 ] with T = 2.0 x 10 9 K described 
by the variable range hopping (VRH) mechanism. 13 ' 
This implies that these systems are in the Anderson in- 
sulator due probably to the strong randomness with the 
localization length £; oc , far longer than the unit cell of 
V(TCNE) X . At that time, we have also to take into ac- 
count the effect of itinerancy on the basis of band picture. 

6. Band Model - Characteristics of Crystal 
Structure 

The feature of the structures of V(TCNE) :r as shown 
in Fig. 2 is that they have planes of V and TCNE with 
rare numbers of pillars of TCNE molecules connecting 
these planes. The electronic state of this 2-d structure 
of V(TCNE)2 will be discussed on the basis of the tight 
binding (TB) approximation, by taking into account only 
d yz , d zx orbitals and LUMO of TCNE as in the case of 
the cluster model. The model Hamiltonian H2d of the 2-d 
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Fig. 6. Two dimensional plane of V(TCNE)2. 

plane of V(TCNE) 2 is written as 

H2d = 12 ( £ i kd (i, 2 .)k^2k + eik*d[i,«)kW»k 

k ^ 

+^2kd (lyz)k (fik + e2k*rf(i ) j /z)k ^4k 



-ei k *dt 



{2,zx)k?^ + £ lkrf( 2 , 2a; )k¥'4k 



+e 2 k*rf (2i!yz)k ^2k + eakdJ 2 ,w a )kV3k 

£ lk d (3, za; )k < ^3k 



^lk* rf (3,^)k^2k 



+zi**d\ Z vz)v ipxv + e 2 kdj 3 ^ )k <^4k + h. 



4 

■ 12 

i=l 



(9) 



where is the energy level of LUMO relative to the d- 
orbitals, and * denotes taking complex conjugation, with 
Qk(* = 1,2) given as 

eik = texp(—ik x a/4: — ikyb/4:), (10) 

e 2 k = -texp(ik x a/A- ikyb/4), (11) 

where a and b are lattice constants as shown in Fig. 6. 
The numbering of the operators of V and TCNE is shown 
in Fig. 6. Here, the terms of d xy orbital are ignored, be- 
cause these localized orbitals are decoupled from H 2 d- 

The matrix of each k of this Hamiltonian (3) is given 
as two matrices of 5 x 5 as follows: 
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One of these eigenvalues of eq.(12) or eq.(13) is zero in- 
dependent of wave vector k and A l ■ Namely, there exist 
the flat bands, since the structure of this matrix is the 
same as the superhoneycomb lattice. 14 '' 

Fig. 7 shows the result of the band structure along 
the symmetry axis in k space and the density of the 
states (DOS) for a choice of parameter of Az,/t = 1.0. At 





DOS(stafeAiinitcell) 



Fig. 8. DOS for the structure of Fig. 2(c). Black and red lines 
represent the cases of t' - Uar /t = 0.0 and t' pi n ar /t = 0.5 where 
t' is the transfer integral between the d yz (or d zx ) orbital of V 
and the LUMO of TONE of the pillar. 



Fig. 7. Band structure of V(TCNE)2. The number in () indicates 
the degeneracy of each orbitals. 



E/t = 0, two degenerate flat bands are found leading to 
the delta function in DOS. The wave function of the flat 
bands are constructed only from d components. Namely, 
of six orbitals of d character, two linear combinations of 
d-orbitals are decoupled from n-d hybridization bands 
such as the NBMOs in the cluster model, which is the 
origin of the flat band. 

In the case of Ad/t > 0, where Ad is the energy level of 
d xy orbital, the Fermi level is located on the flat bands. 
A possible ferromagnetism associated with a flat band 
has been extensively discussed previously. 15 ^ 

In the case of 2-d version of V(TGNE)i,5, a flat band is 
also found on the basis of TB approximation. In this case, 
the occupation of the flat band is incomplete. Therefore, 
the same mechanism as V(TCNE)2 is expected to stabi- 
lize the ferromagnetism. 

In the actual system, however, the flat band is smeared 
to some extent because of three dimensionality and impu- 
rity effects. Nevertheless, the large DOS is still expected 
at the Fermi level. Indeed, from the photoemission spec- 
troscopy experiment, the exsistence of the large DOS of 
d orbital component has been observed near the fermi 
level. 16 ) In the theoretical description, the DOS for the 
structure of Fig. 2(c) has the very same aspects as shown 
in Fig. 8. Here, the effect of pillars of TCNE is taken into 
account properly. Note that there remains huge (techni- 
cally divergent) DOS at the Fermi level, which is consis- 
tent with the experimental result, even for the reasonable 
value of the transfer parameter t' piUar = 0.5i between d yz 
(or d zx ) orbital and LUMO of TCNE of pillar. The origin 
of the peak structure of DOS around the Fermi level can 
be traced back to the flat band in the two dimensional 
sector shown in Fig. 6. Therefore, the ferromagnetism in 
V(TCNE)a; (1.5 < x < 2) can be explained fairly well by 
the mechanism which Kanamori proposed for the ferro- 
magnetism of Ni metal. 17 ' 

In compounds of 1.5 < x < 2, there inevitably exists 



rather strong randomness due to the defect of V atoms. 
This randomness is expected to cause the Anderson lo- 
calization of the quasiparticles at the Fermi level located 
in quasi- flat band. 18 ) Then, the T-dependence of the re- 
sistivity is expected to follow that of the VRH mechanism 
as observed in experiment. 13 -' 

In the case of A d /t < 0, the Fermi level of V(TCNE) 2 
is located on the level of the d xy orbital without disper- 
sion, because the d xy orbital is a localized one. In the 
actual system, however, the d xy orbital forms a narrow 
band, because it weakly hybridizes with p CT orbitals of N. 
Therefore, also in this case, Kanamori mechanism arising 
from hidden flat bands can be considered as the cause of 
ferromagnetism of V(TCNE) X . 

7. Conclusion 

In conclusion, we discussed the origin of the room tem- 
perature ferromagnetism of V(TCNE) X (1.5 < x < 2) on 
the basis of the truncated model of the crystal structure 
newly proposed. We estimated the exchange interaction 
coupling between d-electrons of V atoms by cluster model 
calculations. In parallel, the effect of itinerancy of elec- 
trons on the magnetism was discussed on the band pic- 
ture. It was found that a sharp peak of the density of 
states appears around Fermi level due to the existence of 
the flat bands resulting from the two dimensional plane 
of V(TCNE) 3; (x=1.5 and 2), which gives rise to the fer- 
romagnetism. 

Hence the ferromagnetism of V(TCNE)a; (1.5 <x< 2) 
is theoretically understood both in the localized cluster 
model and the itinerant band picture consistently, which 
will indicate the fact that this ferromagnetism is very 
robust. At the same time, the semiconduting transport 
properties (variable range hopping) in the wide range of x 
is proposed to be due to the exsistence of that band (with 
the large density of state) at the fermi energy affected by 
the randomness. 
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